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We present a theory for the time-resolved optical spectroscopy of high- temperature superconductors at high 
excitation densities with strongly anisotropic electron-phonon coupling. A signature of the strong coupling be- 
tween the out-of-plane, out-of-phase O buckling mode (Big) and electronic states near the antinode is observed 
as a higher-energy peak in the time-resolved optical conductivity and Raman spectra, while no evidence of the 
strong coupling between the in-plane Cu-O breathing mode and nodal electronic states is observed. More in- 
terestingly, it is observed that under appropriate conditions of pump fluence, this signature exhibits a re-entrant 
behavior with time delay, following the fate of the superconducting condensate. 

PACS numbers: 74.72.-h, 71.38.-k, 74.25.Gz, 74.25.nd 



I. INTRODUCTION 



Since its disco veiy in 1986, high-temperature superconduc- 
tivity in cuprates has been a central topic of study in con- 
densed matter physics. It is now widely believed that Cooper 
pair formation is essential for the superconducting conden- 
sate in these systems. However, the nature of the media- 
tor (or glue) responsible for Cooper pairing remains hotly 
debated. Although the interaction of electrons with lattice 
vibrations is not likely to solely account for the essential 
properties of high-Tc superconductors (HTSCs), many probes 
including angle-resolved photoemission,— — inelastic neutron 
scattering^, tunnelingi^ and Raman^ spectroscopies have re- 
vealed that electron-phonon interactions have significant ef- 
fects on various properties. Complementary to the time- 
integrated techniques, different ultrafast pump-probe tech- 



nique: 



,10-15 



have been used to disentangle microscopic inter- 



actions in HTSCs. These techniques aim to study the recom- 
bination of photoexcited quasiparticles and the resulting re- 
covery of the superconducting condensate. In HTSCs, time- 
resolved (TR) angle-resolved photoemission spectroscopy— 
and TR optical reflectivity™ have indicated that the excited 
quasiparticles preferentially couple to a small number of 
phonon subsets before decaying through anharmonic coupling 
to all other lattice vibrations, in support of the notion that se- 
lective optical phonon modes give rise to anisotropy of the 
electron-phonon (el-ph) coupling.- In addition, this anisotropy 
has also been observed in TR electron diffraction, '■ while res- 
onant femtosecond study of both electronic and phononic de- 
grees of freedom suggests strong el-ph coupling^ 

Despite considerable progress in pump-probe experimen- 
tal studies, work on microscopic modeling of the influence of 
el-ph interactions on observables, such as the time-dependent 
optical conductivity or Raman spectra, is very limited. Un- 
derstanding the non-equilibrium dynamics of quantum many- 
body systems has, in fact, posed a theoretical challenge. His- 
torically, theoretical attempts to model the time evolution of 



properties have either used quasi-equilibrium models such 
as T* and /i* models^S to describe non-equilibiium excita- 
tions created by a pump pulse^l or rate equation approaches 
based on the phenomenological Rothwarf- Taylor model- to 
describe the recovery dynamics of the superconducting state. 
Recently, the time evolution of the optical conductivity has 
been studied within a microscopic model that treats the ex- 
citation and relaxation dynamics on the same footing.— All 
these theories are suitable for pump-probe experiments with 
low excitation fluence, where the superconducting condensate 
is merely perturbed but not destroyed. A picture of the dynam- 
ics of quasiparticles and the superconducting condensate in 
the photo-induced phase transition regime.r^— which is im- 
pulsively driven by a high excitation fluence, has as yet been 
beyond reach. 

Here we formulate a theory for the TR optical conductiv- 
ity and TR Raman scattering in HTSCs in the regime of in- 
termediate to high intensity of pump fluence. The theory is 
aimed to address directly the situation where the supercon- 
ducting condensate can be destroyed by pump pulse. It is 
based on an effective temperature model for different sub- 
systems contributing to the response: electrons, hot phonons 
(i.e., out-of-plane out-of-phase O buckUng Big phonons and 
half-breathing in-plane Cu-O bond stretching phonons) that 
are strongly coupled to electrons, and the cold lattice. The 
model phenomenologically includes the effect of the pump 
pulse but addresses in greater depth the electron-hot phonon 
coupling based on a microscopic model Hamiltonian for d- 
wave superconductivity in HTSCs. This microscopic treat- 
ment goes beyond previous effective models for the normal 
stater^i2i allowing us to describe the quasiparticle dynamics 
in both the normal and superconducting states with the same 
approach. Within this unified model, the time evolution of the 
whole set of experimental measurables can be calculated in a 
streamlined way. Our first test of this approach considered the 
Big phonons as the only hot phonon mode in the calculation 
of the TR spectral function for a very high pump fluence j2S In 
the present work, we include the half -breathing phonons as a 



second hot phonon mode in the calculation of the TR optical 
spectroscopy. Importantly, the influence of excitation density 
on quasiparticle dynamics in HTSCs is also investigated. Our 
calculations show that, in the superconducting state, in addi- 
tion to the peak in the optical conductivity and Raman spectra 
due to the Dmde response, there is another peak at higher fre- 
quencies. This high frequency peak disappears when the sys- 
tem evolves into the normal state, but recurs if the supercon- 
ducting condensate is recovered, suggesting the significance 
of the superconducting gap in the TR optical properties. 

The outline of the paper is as follows: In Sec. HI] we lay 
down the effective-temperature model for the d-wave super- 
conductor with electronic coupling to both Big and half- 
breathing stretching phonon modes. The time-dependent ef- 
fective temperatures for the respective subsystems are eval- 
uated depending on the strength of pump fluence. With the 
obtained time-dependence of effective temperatures, the time- 
resolved optical conductivity and Raman spectra and their 
pump-fluence dependence are presented in Sees. Hill and HV] 
Finally, a conduction is given in Sec.lV] 



II. EFFECTIVE TEMPERATURE MODEL 




FIG. 1 : (Color online) Time evolution of effective temperatures of 
electrons Te (red-solid line), Big-phonon mode Tph.i (blue-dotted 
line), half-breathing phonon mode Tph.2 (black-dashed line), and 
cold lattice Ti (magenta-dashed line) for the powers (a) Po = 
0.436 ^W and (b) Po = 0.022 ^W. The inset in panel (a) shows 
a zoomed-in view of Te rising up above Tc, while that in panel (b) 
shows an enlarged view of Te cooling down below Tc- 



Let us consider a two-dimensional superconductor exposed 
to a laser field. The model Hamiltonian can be written as^iZS 






kqzv' cT 



Afield (t), 



(1) 



where c]^^ (bl„) and Ckcr (bi^q) are the creation and annihila- 
tion operators for an electron with momentum k and spin a 
(phonon with momentum q and vibrational mode i/; v = 1,2 
represent the Big and half-breathing modes, respectively), 

Aiyq = ^L -q + bvq, Ck IS the normal-state energy dispersion, 
fi the chemical potential, Ak = (Ao/2)(cosfc2; — cosky) is 
the dx2_y2-wave gap function, Nj^ is the total number of lat- 
tice sites, and gi, the coupling matrix. Following the procedure 
sketched in Ref.|28l we arrive at a four-temperature model: 



9Te 

dr 

dTph,v 
dr 

m 

dr 



= E 



^^X.(Te,Tp„,,) + ^, (2a) 

(2b) 
(2c) 



e L-e 

V 

Ki,{Te, Tph,v) Tph,v — Ti 



C, 



ph.v 



Gr, 



Ci 



ph.u \ -^ph,iy 



T„ 



Ti 



Tp^v 



Here Ki^ is the el-ph coupling kernel, which can be calculated 
from the model Hamiltonian ([U with the equation-of-motion 



approach. It is given by 

47r 
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Nr 



y^g^(ukMk-q ~ 'ykWk-q)^(5(i?k-q - -Bk - ^yci) 



kq 



X n 



l/q 



„(/3ph-/Je)n„q 



./k(l - /k-q)iV0„<, , (3) 



where the Bogoliubov amplitudes are Uk = [(1 + 
i^/E^)/2 Y'^ and v^ = sgn(Ak)[(l - Ck/£^k)/2]i/2 ^j^^ 
ii'k = vCk + Aj being the quasiparticle energy, and the 
Fermi-Dirac and Bose-Einstein distribution functions are 
givenby /k = /(^k) = l/(e'3=^'<4-l), andiVo„ = 7V(f]o) = 
l/(e/3p.no _ 1) ^ith /3,(p,,,,) = l/kBT,(j,h^,y In Eq. ©, the 
specific heat for electrons per unit cell is found to be 









2El + /JeAfc 
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, (4) 



while that for each hot phonon mode is given by 



^ph,v — ~~7~[il^^i^Pph.L>) 



coth 



2 / hllyPpl^jj 

\ 2 



(5) 



in the Einstein mode approximation il^q = Q.^,. Finally, Pf, 
is the power intensity (i.e., power per unit cell) for pump- 
ing electrons and Tfj^ the anharmonic decay time of each hot 
phonon mode. 

Throughout the paper, we use a five-parameter tight- 
binding model2^ to describe the normal-state energy 
dispersion, which is typical to the optimally doped 



BizSrzCaCuaOs+j; (Bi-2212): 
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~2i(cos kx + cos ky) — At' cos k^ cos ky 

-2t"{cos2kr, +cos2ky) 

-4i"'(cos 2kx cos ky + cos kx cos 2/;:^) 



—At"" COS 2/03; cos 2fcy — /.( 



(6) 



where the hopping integrals t = 1, t' = —0.2749, t" = 
0.0872, t'" = 0.0938, <"" = -0.0857, and /i = -0.8772. 
The absolute energy of t is 150 meV. A feature of this disper- 
sion is a flat band with a saddle point at the AI points of the 
Brillouin zone. The temperature-dependence of the d-wave 
gap magnitude is given b}i^ 

AoiT,) = Aaotanh{{Tr /z)^ar{TjT, - 1)} , (7) 

where z — Aoo/(fc_B7c). In our calculations, we set Aqo = 
30 meV, the critical temperature Tc ~ 104 K (from the set- 
ting of Tc = 0.06t for simplicity), the specific heat jump at 
Tc is r = ACe/Ce -- 1.43, and a =^ 2/3. We take flie 
anisotropic el-ph coupling in the form given in Refs. l3ll3ll with 
Oi = 45 meV and g[°^ = 90 meV and ^2 = 70 meV 
and ^2 — 120 meV, T/3,1 — r^.2 = 880 fs, and Cph.i ~ 
Cph.2 = 0.2C;. The pump is represented by a Gaussian pulse 
P = Poe"^'/(^'^'\ with a FWHM (fufl width at half max- 
imum) of 2.35 cr. Hereafter, we set a ~ 4.4 fs, giving a 
FWHM of about 10.34 fs. This value is smaller than the com- 
monly used experimental values of about 30-50 fs but is in- 
deed close to that value of 12 fs used in the recent TR exper- 
iment on YBa2Cu307_5.~ We believe that this variation will 
not affect the qualitative physics, as will be presented below. 
We take the number of k points to be 40 x 40 in the Brillouin 
zone for the temperature evolution and use 256 x 256 for the 
TR optical conductivity and Raman scattering spectral func- 
tion. All calculations are done with the system initially in the 
superconducting state, for T <C Tc. 

Figure [T] shows the time evolution of the effective tem- 
perature for each subsystem for (a) a large pump power in- 
tensity Po = 0.436 /iW and (b) an intermediate value of 
Pa = 0.022 /iW. These values of pump power corresponds 
to 120 fiJ/cxn^ and 6 /.tJ/cm^ of pump fluence in Bi-2212 by 
assuming a 60 nm of optical absorption depth. Starting from 
the initial temperature Tg = Tph^ = Ti = 17 K, the electron 
temperature Te increases rapidly after photoexcitation, and the 
superconductor is driven into the normal state while exhibiting 
a kink structure at Tc (inset of Fig. [TJa)) for both power inten- 
sities. In addition, Tc shows a very narrow peak for the large 
Pq (see Fig. [Ua)). with a broader peak for the intermediate 
Pq (see Fig.[TJb)), due to the fact that the highest temperature 
achieved by electrons is very sensitive to the pump fluence. 
In both cases, the hot phonon subsystems are first heated up 
through their coupling to the photoexcited electrons, and af- 
ter reaching their maximum temperature, they cool down by 
dissipating energy into the cold lattice through anharmonic 
coupling. A noticeable difference between the large and in- 
termediate pump fluences is that in the latter case, the super- 
conducting state recovers more rapidly (i.e., Tc < Tc) in a 
very short time (^ 650 fs), with the kink recurring during the 



cooling stage (inset of Fig. [Ub))- Further energy relaxation 
is then slowed down significantly due to the opening of the 
superconducting gap. 



III. TIME-RESOLVED OPTICAL CONDUCTIVITY 

Within the Kubo formalism, the real part of the TR optical 
conductivity is given by)^ 



ai{LLi) = - 



e^Imn(a;) 



(8) 



where 



I„^H = -|fI:(||)^„(k,.). c9) 



and 



/n(k,a;) = 2i / (ir'TrIm[F*{k,r')A(k, r')]e'"^ (10) 

Here A(k, t') and P(k, r') are the Fourier transform of the 
TR spectral functions i(k, e) and F(k, e) = i(k, e)/(e). In 
the derivation of Eq. (|9]l, we have used the Hilbert transform 



g{'k,iujn) 



-00 «^n - e 



(11) 



with the single-particle spectral function 

i(k, e) = Im[g(k, iw„ ^ e + i5)] , (12) 

n 

where the Green's function g and A are 2 x 2 matrices in the 
Nambu space. Since this spectral function as calculated with 
the method of Ref. |28| is a function of the effective electronic 
temperature, which is time dependent (see the discussion in 
Sec.|II]i, it is time resolved. Therefore, the optical conductivity 
and the Raman spectra as discussed in the next section are also 
time dependent. 

Figure[2]shows the time evolution of the real part of the op- 
tical conductivity 0-1(0;) at several selected time delays. From 
Fig. 121 one can see that at all time delays, the optical con- 
ductivity shows the well-known Drude peak at a; = 0, due 
to the nodal quasiparticles for the d-wave gap symmetry. In 
addition, for r = —440 fs and —11 fs, at which the mate- 
rial is superconducting and Ao « 30 meV, we observe that 
(Ji{lu) exhibits a broad peak at about uj — 2Ao + ^ph,i- 
(The specific location may be affected by several factors, al- 
though 2Ao + rip/i 1 plays a substantial role.) Our obser- 
vation is consistent with earlier study of optical conductiv- 
ity in the thermal equilibrium state of HTSCs.^^^ In contrast, 
no peak at u; = ^ph.2 is observed (a signature of the cou- 
pling between electrons and the half-breathing mode along 
the nodal directions). This is due to the fact that the coupling 
between electrons and the Big phonon mode is the strongest 
at the M points, at which the van Hove singularity is also 
located. This is further verified by the observation that no 
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FIG. 2: (Color online) Real part of the time-resolved optical conduc- 
tivity at several different time delays for (a) Pq = 0.436 /iW and (b) 
0.022 fiW. The number in parentheses is Te. 
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IV. TIME-RESOLVED RAMAN SCATTERING SPECTRUM 

The time-resolved Raman scattering intensity is calculated 
via a simple relation^ from the imaginary part of the Raman 
response function. In the bare vertex approximation.r^i^ it is 
found as: 



2n^ 
Imxiinm -^UJ + i5) = -^rY. 7k^x(k, w) , (13) 



where 



/x(k,w) = 2i f dr' Trlm[f3i^*(k,r')f3i(k,r')]e*"^'. 

(14) 
Here 7^ is the nonresonant bare Raman vertex given by 
7k = 7k "Ta with fa being the Pauli matrix and 7k = 



Z^a,/: 



.S_d^ 



./,S 



are the polarization unit vectors 



120 

to (meV) 



of the incident and scattered photons, and ^k the electronic 
normal-state dispersion of the conduction band. 

Figure [3] shows the time evolution of the Raman scattering 
spectrum. When the electron-hot phonon coupling is switched 
off, the Raman spectrum rises with cj and has a large peak at 
twice the gap, 2Ao, at the initial temperature (red solid line 
in Fig. [3]). A small shoulder in the curve around 90 meV 
arises due to the van Hove singularity. In the presence of the 
electron-hot phonon coupling, the superconducting gap func- 
tion is renormalized, which shifts the original 2Ao-peak to 
lower frequencies. Simultaneously, another peak develops at 
2Ao + ^ph.i but no peak develops at rip/1,2, for the same rea- 
son for the optical conductivity (Fig. |2]i. After photoexcita- 
tion, this double-peak structure evolves into a very broad peak 
as the system enters the normal state. For large Pq, this broad 
peak remains for a few picoseconds. However, for intermedi- 
ate Pq, once the superconducting state recovers, the double- 
peak structure appears again. This result is fully consistent 
with our calculations of the TR optical conductivity described 
above. 



FIG. 3: (Color online) Time-resolved Raman scattering spectrum 
at several different time delays for (a) Po ~ 0.436 /iW and (b) 
0.022 /.tW. The number in parentheses is Te. 



such peak appears when the el-ph coupling is turned off in 
the initial superconducting state Te = 1 7 K (red solid line in 
Fig.|2]). After photoexcitation, the superconducting gap, with 
its time-dependence encoded in the effective electron temper- 
ature Ao(Te), is decreased and the high-frequency peak shifts 
toward lower frequencies, merging into the zero-frequency 
Drude peak in the normal state. For large Pq, the Drude 
peak remains for the whole time-delay range being simulated. 
However, for intermediate Pq, once the superconducting con- 
densate recovers, the (2Ao(Te) + l^p/j i)-peak recurs — a re- 
entrant behavior The absence of the peak at nph,2 is robust 
with time delay. 



V. CONCLUSION 

We have presented a theory for the time-resolved optical 
conductivity and Raman spectrum, based on the TR spectral 
function that we have recently formulated for HTSCs. Our 
calculations show that the signature of the electron-i?ig mode 
coupling in the TR optical conductivity and Raman spectrum 
is more pronounced than the consequence of the coupling be- 
tween electrons and the half breathing mode. This is the result 
of a concurrence of anisotropy of the el-ph coupling, band 
structure, and d-wave energy gap in HTSCs. Even more in- 
terestingly, this signature also shows a re-entrant behavior in 
concurrence with the superconducting condensate, which can 
be controlled by the pump fluence. The observation of the 
broad peak in the TR Raman spectra and their re-entrant be- 
havior provides a direct evidence of the el-ph coupling. 
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